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Abstract: 

Various methods are discussed for obtaining the CKM angle 7 
through the interference of the charged i?-meson decay channels 



B- K-D^ and 5" ^ R- where the and decay to 
common final states. It is found that choosing final states which 
are not CP eigenstates can lead to large direct CP violation which 
can give significant bounds on 7 without any theoretical assump- 
tions. If two or more modes are studied, 7 may be extracted with 
a precision on the order of ±15° given ~ 10^ 5-mesons. We also 
discuss the case of three body decays of the D where additional 
information may be obtained from the distribution of the D decay 
products and consider the impact of DD oscillations. 
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1 Introduction 



The only manifestations of CP violation observed to the present time are 
those in the neutral kaon system. The advent of several machines capable 
of producing large number of B mesons makes it likely that at least some 
examples of CP violation in B physics will soon be seen. 

B factories at SLAC and KEK have been specifically constructed to ob- 
serve CP violation in the time dependent oscillation of B^B when B^ decays 
to a CP eigenstate such as ipKg or 27r. In such experiments one has the ad- 
vantage that the CP violating phase may be in a few cases cleanly extracted. 
In contrast direct CP violation which can occur in B^ decay always appears 
in combination with a strong phase which cannot be easily determined since 
it has its origin in strong interaction physics. In order to put the extraction 
of the CP odd phase on the same footing as the clean oscillation experiments 
it is thus important that some way of eliminating the uncertainties due to 
the strong phase be found. 

Direct CP violation in B^ decay may prove to be an important component 
of the future B physics program because it offers signals sensitive to the 
angle 7 of the unitarity triangle |l|] . Here we focus on effects which originate 
from the interference of B^ K^D^ with B^ K~D^. In contrast, 
experiments involving B^ may only extract a (from B^ 2tt for instance) 
and (3 (from 5° — >■ ipKs)- The unitarity triangle which follows from the 
unitarity of the CKM matrix is a key prediction of the Standard model. 
Independent measurements to over determine the triangle by measuring each 
of its sides and angles therefore provide a non-trivial test of the standard 
model. Furthermore, the Standard Model predicts that a number of different 
measurements of CP violation will depend on the same phase 7 (e.g. direct 
CP violation in i? ^ Ktt and oscillations in Bg). Comparing the results of 
these experiments will thus be a sensitive test for new physics. 

Since the modes we consider are direct and not time dependent, they may 
be observed in any experimental setting where large numbers of B mesons 
are produced. Aside from the SLAC and KEK asymmetric 5-factories, these 
include CLEO and hadronic B experiments such as BTeV, CDF, DO and 
LHCB or a high luminosity Z factory. 

A method for putting the determination of 7 on the same footing as the 
oscillation experiments by using direct CP violation was first suggested in P| 
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where the decay B ^ K is considered followed by the decay of to 
a CP eigenstate (see also [0]). In this case the fact that the final state is a 
CP eigenstate means that it will interfere with the channel B~ -h> K~D^ . In 
Section 2, for self containment, we will review the important features of this 
method and discuss some problems which arise in its implementation. 

These problems can be resolved through a more general method using 
non-CP eigenstates suggested in 0, |^ and refined in ||^ which not only enable 
a clean extraction of 7 but in addition have the attractive feature that they 
can give rise to large CP asymmetries. This is outlined and expanded upon 
in Section 3 and 4. In Section 5 we estimate some relevant branching ratios 
and obtain a rough estimate of the attainable accuracy in extracting 7. 

In Section 6 we discuss the somewhat more general case where the dis- 
tributions of three body decays of are considered. Here we consider two 
different methods for obtaining 7. First, one can fit the distributions to a 
resonant channel model where each of the channels can be considered a quasi 
two body mode. The phases between the channels will thus give additional 
constraints on the fit. Second, if each point on the Dalitz plot is considered 
a separate mode, in some cases the accumulated inequality bounds on sin^ 7 
can provide a determination of sin^ 7. In section 7 we give our summary and 
conclusions. 

Throughout this paper we will assume that the effects of DD mixing 
are negligible. In an appendix, we discuss the impact of DD mixing on 
this method for determining 7. Specifically, if one uses the time interval 
between the parent B~ decay and the subsequent D decay, one can eliminate 
the effects of mixing. One can also obtain information about 7 in time 
independent studies. 



2 Using CP Eigenstates Decays of 

Gronau London and Wyler's strategy (GLW) to extract 7 from CP violation 
in the decay B^ — > K^D^ followed by D^ — > CP eigenstate is to separately 
determine the branching ratios |]^: 



(a) Br{B- R- D 

(b) Br{B- K-D 
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(c) Br{B- K-Dl) or Br{B- R- D^_) 

together with their conjugates where denote the CP eigenstate = 
(Do±^o)/v^. 

Once a, 6, c and the corresponding quantities for the conjugate modes are 
known then one can separate out the interference effects and thus determine 
cos(Ca' + 7) and cos{C,k — 7) simultaneously where C,k is the CP conserving 
strong phase difference between B^ K^D^ and B~ — > K~ D while 7 is 
the CP violating weak phase difference. From these two cosines, the values of 
the actual phase angles C,k and 7 can clearly be determined, up to a eight-fold 
discrete ambiguity as will be discussed in more detail in section ^ 

The rate for (c) is experimentally observed through the decay to a CP 
eigenstate such as K~^K~ , Ks4> and -7r"'"7r~ and presents no problem in prin- 
ciple. Likewise the decay (a) should be readily measurable through ei- 
ther leptonic or hadronic modes of the D^. As it stands, however, this 
method has a very serious problem in measuring the branching ratio of (b), 
Br{B~ — > K~D^). The detector must distinguish from to determine 
the decay rate to this rare mode. The background decay Br{B~ — > K~D^) 
(i.e. (a)) is expected to be larger by about two orders of magnitude. 

There are only two ways that could possibly be used to tag the flavor of 

the 

(a) through semi-leptonic decays 

(b) through hadronic decays 

We will consider each of these individually and show that neither is fully 
satisfactory. 

The semi-leptonic tag, i.e. the quark level process c -^"sz/^, has the 
problem that there is an overwhelming background from the direct semilep- 
tonic decay of the B meson. The signal here is the sequential decay: 

B- R-Jf ; Tf rVeXj (1) 

while the background is from the direct semi-leptonic decay of the b quark 
followed by hadronic decays of the charm quark: 

B' i'VeXc ; sud (2) 
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Both give rise to the same sign lepton and while there are several features 
distinguishing the signal from the background it represents a serious problem 
as the branching ratio for the signal is expected to be O{10^^) whereas that 
for the background is expected to be C(10~^). The signal/background ratio 
is thus daunt ingly small. 

Let us note some of the characteristics of the signal which distinguish it 
from the background. First of all, in the B~ rest frame the K~ tends to be 
monochromatic. Also the signal tends to yield a pair of kaons K~K^ (or 
K~K^) whereas the background dominantly has one kaon, K~ or K . In the 
signal the semi-leptonic decays of the originate from a tertiary vertex in 
sharp contrast to the case of the background. It is difficult to say whether 
these distinctions are enough to separate the signal. Detector specific studies 
are required to give reliable answers but the size of the signal/background 
ratio does not give us much ground for optimism. 

While the use of such a leptonic tag is likely to be impractical, the 
hadronic decay has a more fundamental problem. In this case one would 
detect the decay B~ K~D by observing a final state where the D decays 
through a Cabibbo-allowed (CBA) process, for instance K^tt". Un- 

fortunately, the doubly- Cabibbo suppressed (DCS) decay of D° (e.g. — >• 
K^n^) will also lead to the same final state at a branching ratio two orders 
of magnitude smaller. 

On the other hand, as was pointed out in [^, the initial decay B~ — > 
K^D is color-suppressed (CLS) while the decay B~ — ^ K~D^ is color- 
allowed (CLA). Thus difference in the production rate tends to offset the 
difference in the decay rates for the two processes and since the final states 
are really indistinguishable, they will interfere quantum mechanically. 

As a specific example consider the possible tag K^tt~ . The signal 

from B^ —>■ D K~ and the background from B^ D^K^ will be given by 
the sequences: 



A) B- K-D ■ D K+7V 

B) B- K-D^ ; D° ^ fsT+TT 

These two decay chains have the qualitative form. 



(3) 
(4) 



A: 



CLS ® CBA 



(5) 
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B : CLA ® DCS 



(6) 



Numerically the ratio of these two amplitudes is appreciably close to 1 since 
it is expected to be given by: 



M{B- K-D K+TT-]) 



Oi 
^2 



B{Tf K+TT- 



(7) 



Here ai and 02 control the relative sizes of the CLA and the CLS amplitudes. 
Experimentally [0 the indications are that 



02 
ai 



0.26 ±0.07 ±0.05 



(8) 



which roughly agrees with the simple color counting value of 1/3. 

Let us consider B[D^ K^-n^) / B{D^ K^-n^) which is of order A^. 
We can formulate an estimate taking into account SU(3) breaking effects in 
the form factors and the decay constants. For the present purpose we use 
the experimental result §]: 



Br{D° 



TT 



Br(n^ K+n- 
as well as the ratio of CKM elements 



.0072 ± .0025 



ub 



cb 



.08± .02. 



(9) 



(10) 



Then with A = 0.23 and the central values from (Eqns. H-10) we therefore 
get: 



M{B- 



TT 



1) 



(11) 



The two amplitudes are roughly comparable and we cannot tell whether the 
charmed particle was a or a D . B~ — > K~D with a D decaying 
hadronically will give rise to a final state which is indistinguishable from the 
corresponding decay of the in B^ — K~D^. The two amplitudes will 
thus be subject to large quantum mechanical interference effects. The use 
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of a hadronic tag for determining B{B^ K^D ) for the GLW method 
appears therefore to be ruled out. 

Despite these difficulties, one need not discard the GLW approach. The 
only input which is lacking is the branching ratio K~D^ . It may be 

possible to theoretically estimate this quantity which will allow the GLW 
program to go forward. 

Here we consider an alternative approach where we take advantage of 
these interference effects to enhance CP violation and ultimately provide 
another way for a clean (i.e. no penguin pollution) way of extracting 7. As 
discussed in ||^ the fact that these two amplitudes have large interference 
effects implies that there will be large CP violating asymmetries in such 
combined decay channels which in turn gives us a handle on measuring 7. 
Thus, it is instructive to consider why CP violation will be enhanced in this 
case as compared to GLW. 

In the eigenstate case, the size of the expected CP asymmetry is controlled 
only by the ratio of the amplitudes for B~ — > K~D^ versus B~ K~ D . 
Following the above estimate and taking into account the appropriate CKM 
factors: 

M{B- K-lf) 
M{B- K-D^) 

This means that the maximum possible size of the CP asymmetry is expected 
to be 0(10%). In contrast eq. (pT]) imphes that the two interfering amplitudes 
have roughly the same magnitude and so the interference effects, and in 
particular CP violation, will be near maximal. 

Of course the search for large direct CP violating signals is interesting in 
its own right but the real goal is to extract the angle 7 from experimental 
results cleanly, that is without any reference to a model for hadronization. 
While the original idea of 0, though sound in principle is probably not 
practical, the basic concept may be used if one observes two or more dis- 
tinct hadronic states that are common decay products of and D (as all 
hadronic final states of are). With this information, one can reconstruct 
7 cleanly. We now consider a number of methods based on this idea. 



02 



V,. 



ub 



VcbV: 



0.1 



(12) 
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3 Non-CP Eigenstates Decays of 

Using this idea, let us now consider the case where the two channels 



B 



R-D 



,0 



B 




(13) 



interfere because both and D decay to some common final state X. In 
the GLW method the specific case where X is a CP eigenstate such as KgH^ 
was chosen while we will focus on the instance, considered in 0, where X is 
not a CP eigenstate. In particular, following the logic of the previous section, 
the case where — X is a DCS decay and — > X is a CBA decay is of 
particular interest, for instance X = K^n^ . 

In order to formulate an expression for the rates, let us define the following 
branching ratios: 



a{k) = Br{B- fc" 




a{k) - 


= Br{B+ - 


-^k^W) 


b{k) = Br{B- r 




b{k) -- 


= Br{B+ - 


^ k+D°) 


c(X) = Br{D^ - 


.X) 


c(X) 


= Br(D° - 


-.X) 


c(X) = Br{D^ - 


.X) 


c(X) 


= Br(D^ - 


-.X) 



d{k,X) = Br{B- k-[X]) d{k,X) = Br{B+ A;+[X]) 



Here k^ represents either or K*^ (or indeed one may consider any 
other kaonic resonance or system of strangeness=— 1 and well defined CP) 
and [X] is the common and decay channel observed. Thus the com- 
bined rates d{k, X) and d{k, X) include the effects of the interference of the 
two channels. 

In the standard model, it is expected that a{k) = a{k), h{k) = h{k) and 
c(X) = c(X) all of which we will assume from here on. In general, however 
if 7 7^ one expects d{k,X) ^ d{k, X) (CP violation) and indeed the value 
of the quantities d, d may be expressed in terms of a, b and c as: 



(14) 



d{k,X) 



a{k)c{X) + b{k)c{X) 
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+2^ a{k)h{k)c{X)c{X) cos(a + 5x + l) 
d{k,X) = a{k)c{X) + b{k)c(X) 



+2^a{k)b{k)c{X)c{X) cos(a + - 7) (15) 

where Cfc is the strong phase difference between —>■ k^D^ and B^ — >■ 
k~D^; 6x is the strong phase difference between D X and D ^ X 
and 7 is the CP violating weak phase difference between B^ k^D^ and 
B- k-D^. 

In the standard model 7 is given directly from the CKM elements 



l = arg{-Vu,V:,V,,V:,). (16) 

Existing data does not constrain 7 very much giving an allowed range in the 
SM at 95% c.l. [I of 36° < 7 < 97° corresponding to 0.35 < sin^ 7. 

The strong phases (k and 6x that result from hadronic final state inter- 
actions cannot be reliably calculated with any known method and must be 
determined experimentally. Here we will take the approach that information 
about 6x and (k are extracted from the data along with b{k) and 7. 

The above may be made somewhat more general by considering the class 
of modes B^ k~D where D is an excited D meson. Let us suppose that 
D subsequently decays D —>■ D + where is a single particle and the D 
then decays into a CP non-eigenstate mode of the type we have considered 
above. For instance D may be a with k- = K- and we may use either 
of the following decay chains: 



B- K-{D* ^ n^[D ^ X]) 

B K-{D* ^-f[D ^ X]) (17) 

Here — X is a DCS mode. In these sort of examples the analysis would 
be essentially the same as considered for the ground state of the D. The only 
constraint on this generalization is that one of {/c, D} should have spin so 
that there is only one partial wave otherwise multiple partial waves would 
have to be separated and the analysis for the extraction of 7 would be more 
complicated and has been considered in [|T0|. 
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4 Methods for Extraction of 7 



Let us now turn our attention to the extraction of the weak phase 7 and 
incidentally also the total strong phase ^{k,X) = Cfc + ^x- We will consider 
two scenarios under which such a reconstruction is possible assuming that 
the values of c{X) and c{X) have been determined in advance and that the 
rate of D^D^ mixing is negligible. In the appendix we show how the possible 
effects of mixing may be removed from the data. 

Case 1: For one particular mode, {k,X), d{k,X) and d{k,X) are known 
and in addition a{k) and b{k) are known. 

This is the simple generalization of the GLW method to the case where 
X is not a CP eigenstate. We thus assume that the branching ratio b{k) may 
be obtained from tagging the D with the semileptonic decays as discussed 
in section |^ or as the result of some theoretical estimate |Tl|. If b{k) can be 
obtained in some way, this method is less involved than the other method 
we will discuss below. 

As in the GLW method we are required to extract the interference terms 
in d and d by solving the two equations (^) which we can rewrite: 



d{k,X) = 
d{k,X) = 

where R{k, X) is 



a{k)ciX) + bik)ciX)\ [1 + R{k, X) cos(a + 5x + 7)] 
a{k)c{X) + b{k)c(X) \ [1 + R{k, X) cos(Cfc + Sx - 7)] 



(18) 



Defining 



R{k,X) 



2^Ja{k)b{k)c{X)c{X) 
a{k)c{X) + b{k)c(X) 



(19) 



Ai 
A2 



d{k,X) 



R{k,X) [a{k)c{X) + b{k)c(X) 

1 r d{k,X) 
Rik,X) [a{k)c{X) + b{k)c(X) 



=- - 1 



(20) 
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the solution for 7 and ^{k, X) = (k + is: 

7 = -(cr cos~^ Ai — r cos~^ A2) + nvr; 
^{k,X) = i(acos"^Ai + rcos"^A2) + n7r. (21) 

In the above a, r G {±1} and n G {0, 1} expressing the fact that there is a 
eight fold ambiguity [jl2| since the sign in front of each of the cos~^ functions is 
undetermined and we can add tt simultaneously to C, and 7 without changing 
the results. Specifically, the eight solutions giving results identical to a given 
(^,7) are {(^,7), (-'C, -7),(7, 0> (-7,-0^(7^ + ^,7^ + 7), (7r-^,7r-7), 
(tt + 7, vr + ^), (tt — 7, 77 — ,^)}. To resolve the ambiguities between the strong 
phase and the weak phase we can use this method on two or more modes 
since 7 will be the same for each mode while C,{k,X) should be different. 
This is a simple generalization of the GLW method discussed in section in 
that case, c{X) = c{X) while 6x = nn. 

More generally, if we do consider two decay modes with different strong 
phases (mod vr), we can dispense with the need to know the value of b{k) as 
follows: 

Case 2: For two distinct modes {(fc, Xi), (fc, X2)} the quantities d{k,Xi) 
and d{k,Xi) are known. In addition a{k) is known but not b{k). 

In this case we assume that one cannot easily obtain the branching frac- 
tion b{k) using the semileptonic tag or by any other means. However, if d and 
d are known for two different modes, we can solve for the missing information 
(7 and b{k) and ^) up to a discrete set of ambiguities. 

Specifically, for decay modes Xi and X2 we assume that we have 
measured the quantities {a{k), c(Xi), 0(^2), c(Xi), c(X2)} as well as d{k, Xi) 
and d{k,Xi). There are therefore four unknowns that must be solved for: 
{b{k), ^1, ^2, 7}- To do this we use the four equations: 

d{K,Xi) = a{K)c{Xi) + b{K)c(X^) + 2^ a{K)b{K)c{Xi)c{Xi) cos(ei + 7) 
d(X,Xi) = a{K)c{X^) + b{K)c{X^) + 2^ja{K)b{K)c{Xi)c{Xi) cos(ei - 7) 
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(22) 

d{K,X2) = a{K)c{X2) + h{K)c{X2) + 2^ a{K)h{K)c{X2)c{X2) cosfe + 7) 

d{K,X2) = a{K)c{X2) + h{K)c{X2) + 2v/a(ir)6(ir)c(X2)c(X2) cosfe - 7) 

(23) 

To solve these, let us define the quantities. 

_ Q!(fc,X,)+d(fc,X,) _ . . 
^» - 2a(fc)c(X0 

2(zi - Z2)ui + (1 - p)u^; 
= sin^ 7; (24) 

where Zi, e and p are known directly from experiment and Ui and Q = 
sin^ 7 must be solved for. 

The equation which ui must satisfy is easily derived: 

Aui6e = (e - 6){yj6 - {z^ - ui)h) (25) 

Since 6 is second order in ui this equation is in general a quartic equation 
which may have up to 4 real roots. For each real root {ui)k (where k = 
1, . . . , 4 indexes the solutions of eqn. (|25|)), sin^ 7 is given by: 



Ui 



_ b{k)c{Xi) . 
~ a(fc)c(X,)' 

c(Xi)c(X2) ^ 
c(Xi)c(X2) 



d(k,Xi)-d(k,Xi) . 
2a(fc)c(Xi) ' 



M2 ■ 

Ml ' 



_ zl/p - 

yl/p\ Q 



sin'^ = Q = (26) 

e — 

where 6 is given in terms of Ui by Eq. (|2^). Each root leads to a four fold 
ambiguity in the determination of 7; taking up to four roots together there 
is therefore < 16 fold ambiguity in the determination of 7. To reduce this 
ambiguity (especially if it should turn out that all 16 possibilities manifest), 
it is therefore helpful if observations are made of at least 3 modes of 
decay (for a given k) in which case only an overall four fold ambiguity in 
7 remains. Since this method determines Q = sin^7 it, therefore, cannot 
distinguish between the solutions {±7, vr ± 7}. 

For each solution the corresponding total strong phase difference C,i is 
then determined without further ambiguity by: 
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Another way to resolve the discrete ambiguity is to determine indepen- 
dently the phase difference: 



Ae = ^2-^i = Sx,-Sx, (28) 
from the study of D'^ decays |T3]. It is related to the value of 7 and Ui 



through equation (pTl) and in general only two values of 7 and Ui will give 
the correct value of A,^, the remaining two fold ambiguity being between 7 
and 7 + TT. 

To qualitatively understand the solutions of these equations, it is useful 
to consider a plot of 7 versus b{k). First, let us assume that we have perfect 
experimental information. For a given decay mode Xi where we know {a{k), 
c{Xi), c(Xj), d{k,Xi), d{k,Xi)} while j^,, 7, b{k)} remain unknown. The 
two equations (p2D for the mode Xi give a locus of points in the 7 — b{k), or 
equivalently the — Ui plane when C,i is eliminated. Let us now consider the 
properties of this curve. 

Inspection of these equations shows that they are left unchanged under 
the transformations: 



(7,e.) ^ (-7,-0 

(7,e^) ^ (7r + 7,7r + (29) 

Thus it follows that the curve is periodic with respect to 7 — ^ 7 + vr and that 
the curve is also symmetric with respect to 7 — > vr — 7. The curve in the 
range < 7 < 7r/2 can therefore be reflected through 7 = and 7 = 7r/2 to 
get the entire curve. 

For 7 = nn, the two cosines in eq. (^) are the same so if there is any 
CP violation the equations are inconsistent. This is obvious from the physics 
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since 7 is the only CP violating parameter; thus 7 = nvr implies there is 
no CP violation. Conversely, if a finite amount of CP violation is observed, 
some bound can be placed on 7. In particular, a lower bound Qmin can be 
placed on Q: 



and a'{k,Xi) is the CP asymmetry defined by: 



^ ' dik,Xi) +d{k,Xi) 1 + Zi ^ ^ 

We also define jmin to be the angle in the first quadrant such that Qmin = 
sin^ 7min- Since 7mm represents the extreme left edge of the figure eight curve, 
there is a unique value of Ui which gives Q = Qmin we will denote this by -Uj 
which is given by: 



iti = Zi + 2{1 - Qmin) (32) 

Likewise, there can be no CP violation if m — > or u ^ cxd and so the 
observation of CP violation implies an upper and lower bound on the value 
of u and hence b: 



Umax = {I + y Zi - \yi\ + 1)'^ 

Umin = (1 - \/zi + \yi\ + 1)^ (33) 

which leads to the bounds bmin < K^) — bmax where hmin = o,{k)c{X)uminl c{X) 
and hmax = a{k)c{X)umax/ c{X). 

Eliminating from the two equations (|2^ ) gives a quadratic equation for 
Ui for a given value of 7 which defines the set of solutions in the 'y — Ui plane: 

Qm^ - 2{zi + 2(1 - Q))Qui + {Qz} + (1 - Q)y}) = 0. (34) 
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This is quadratic in uf, therefore, there are zero, one or two solutions for 
Ui at any given value of 7. In particular, if Q < Qmin there are no solutions 
while for Q = Qmin there must be exactly one solution (i.e. u = u). It is 
also the case that when Q = 1, for instance, when 7 = there is again 
exactly one solution Ui = Zi. This follows from the fact that in the sum 
d + d the interference term vanishes and a value of Ui may thus be obtained. 
Taking the 7 axis horizontal, the curve in the range < 7 < vr therefore has 
the topology of a lazy eight centered about the vertical line 7 = 7r/2 which 
crosses itself at {'K/2,Zi). 

In Fig. 1 we show a plot of Ui versus 7 for Zi = 1.5 and i/i = 0, 1, 
2. The bounds given by Eqns. ( ^^3]) are indicated with the rectangular 
boxes. Clearly the greater the value of yi (and the greater the amount of CP 
violation), the greater the bounds which may be placed on Ui and 7 through 
these inequalities by considering just one mode. Thus single modes which 
have a high degree of CP violation are quite desirable since they lead to the 
most restrictive bounds in the 7 — ^4 plot. Indeed, as we have argued above, 
larger CP violation is more likely to arise in the case of non-CP eigenstate 
final states (such a.s D ^ K^n^) as compared to CP-eigenstate modes. 

Let us now turn our attention to the case when we have two modes 
present. If we had perfect data concerning each mode, we would generate 
two of these lazy-eight curves in the 7 — b{k) plane. In general, these curves 
will therefore intersect in as many as four points in the range < 7 < 7r/2 
which correspond to the solution of eq. (pSj) . If the data was inconsistent, 
then the curves would miss each other corresponding to a situation where 
eq. ( p5| ) has no real solutions for 7. In order to understand how this might 
play out, let us now consider a scenario for the not yet observed branching 
ratios. 

5 Attainable Accuracy 

5.1 Estimates of Branching Ratios 

For the purpose of illustrating the ability to extract 7 by combining mea- 
surement of several modes, we need to make a simple estimate of the branch- 
ing ratios of DCS decay modes of D which have not yet been measured. 
We will proceed by relating these modes to the well measured Cabibbo- 
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allowed ones and then use factorization to break down D M1M2 to 
{M2\Jcd\D){Mi\Jus\D). The SU(3) breaking of the decay constants is used 
in keeping track of the piece {Mi\Jus\0) . Single-pole dominance allows one 
to keep track of the SU(3) breaking in {Mi\Jcd\D) . In the final estimated 
rates, we also factor in the small difference in phase space. 

For a concrete example let us consider the DCS mode — > K~^7r~ 
which we relate in this procedure to CBA counterpart — »• K'n'^. Using 
this reasoning the amplitude will be proportional to: 

A{D'^K^n-)ocX'fKiml-ml)(l-^] ' (35) 

Thus 



In this instance there is no phase space correction. Using A = 0.22, fx = 160 
MeV, = 132 MeV and other masses we thus get. 



= 5.3 X 10"=^ (37) 
Therefore BR{D^ n+R-) = 3.83 x 10"^ g gives: 

BR{D^ ^ K+71-) = (2.0 ± .7) X 10"^ (38) 
which should be compared to the measured value 

BR{D^ ^ K+TT') = {2.8 ±0.9) xlQ-^ (39) 

We will normalize our predictions for all the DCS modes to the measured 
central value. 

For the other DCS modes of interest to us we can now use relations 
amongst DCS modes and their CBA counterpart. Thus, e.g., for — »• K~^p~ 
we should have: 
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A{D^ K+p-) A{D^ 71+K*- 



A{D^ K+7T-) A{D^ TX+K-) 
Similar scaling relations for the other two modes are: 

A{D^^K+a^) _ ^7r+i^f(1270)) 

A{D^ K+TT-) ~ A{DO -> n+K~) 

A{D^^K*+7i-) _ A{D^^p+K-) 

A{D'^ K+TT-) ~ ^ TT+K-) 



(40) 



(41) 



The resulting branching ratio for all the DCS modes of interest are given in 
Table |. 

To complete our sample calculation, we need to estimate the expected 
magnitudes of {a,b} 

Starting with a{k), we will extrapolate from the observed branching frac- 
tion for related B decays Q: 

Br{B-^7i-D°) = (5.3±0.5) X IQ-^ 
Br{B-^p-D°) = (13.4 ± 1.8) X 10^1 (42) 

Multiplying this by sin^ 9c one obtains the estimates for a{k): 

a{K) ^ 2.6 X IQ-^ 

a(K*) ^ 6.6 X IQ-^ (43) 

The estimation of b{k) is more uncertain since it is a color suppressed 
process. Thus, to estimate the branching ratio for B~ — > k~D^ we use the 
fact that the quark level diagram for this process is color suppressed with 
respect to B^ — > k^D^ and take the color suppression to be simply a factor 
of l/Nc- Folding in all the appropriate CKM elements an estimate for the 
branching ratio may be obtained from the previous estimate of a{k) (taking 
iVc = 3): 

1 2 

^ .015 (44) 







Vub\ 




a{k) 






\Vus\ 
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and so for the two specific cases: 



b{K) ^ 4.0 X 10"^ 

b{K*) ^ 10.0 X 10"^ (45) 
5.2 CP Violation in Single Modes 

Let us now perform some sample calculations using the above estimated 
branching ratios in order to illustrate what precision might be obtained by 
this method. 

First, we will consider what can be learned through the use of data in a 
single mode. As discussed above, this will not allow one to obtain an exact 
value of 7 but can give a bound on 7 if large CP violation is present. We 
will therefore consider how many B decays are needed to have a measurable 
signal of CP violation and what bounds on 7 would follow. We will next 
consider the precision for extracting 7 by combining the data from several 
modes. 

To illustrate this let us consider the specific case of CP violation when 
X = K^Ti~ . The observed branching ratios for the D decays are: 



c(/r+7r-) = (2.9 ± 1.4) X 10"^ 

c(ir-7r+) = (3.83 ± 0.12) X 10"^ (46) 
Where the partial partial rate asymmetry is given by: 



which from equation ([ISD can be written: 

'(h X\ = ^(fe,^)sin(Cfc + (^x)sin7 

l + i?(A;,X)cos(a + 5x)cos7 ^ ' 

Let Nb be the total number of -B^. In a given mode, let Ei be the 
acceptance times efficiency of a detector and let us denote: 
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Nb = EiNB 



(49) 



Thus if N^{k,X) is the number of charged i?^ required to observe the 
partial rate asymmetry at a 3a level with an ideal detector, it is related to 
the actual number, A^|^(A;, X), of by EiNl^{k, X) = N%^{k, X). In terms 
of the CP asymmetry, is given by: 
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^ [a'{k,X)Y[d{k,X) + d{k,X)] 

In order to obtain a large value of a' it is clearly necessary to have a 
large value of R{k,X) since \a'\ < R{k,X). As defined, < R{k,X) < 1 
where R{k,X) is maximized when a{k)c{X) ^ b{k)c{X). It should be clear 
that this will happen if the two channels have roughly the same amplitude. 
In particular, using the numbers in the estimates above, R{K^ , K^n^) fa 
R(K*~ , K'^7r~) ^ 0.94. On the other hand if we had considered the case 
where decays in a CBA mode, then we would have obtained a much 
smaller value, R{K~ , K~7r~^) w R{K*~ , K'tt'^) .02 and so CP violation 
would be small. 

To completely specify a' and N^'' , of course we also need to know cos ^ cos 7 
and sin ^ sin 7. Since these are totally unknown, to get a rough idea of the 
experimental requirements we will take the sample values cos(^)cos7 = 
and sin(,^) sin 7 = 1/2. 

In this example we obtain 



a'{K-,K+7r~) a'{K*-,K+7r-) ^ 0.47 

Nl^'iK-^K^TT-) ^ 17.6x10^ 
Nl''{K*-,K+7r-) f« 7.0x10^ (51) 

Using the estimated branching ratios in Table 1, we can perform a similar 
estimate for some of the other possible modes: 



Nl''{K-,K+p-) ^ 6.3x10^ 
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^ 2.5 X 10^ 






«r) 


^ 9.3 X 10^ 






«r) 


^ 3.7 X 10^ 








^ 13.6 X 10 




K*+ 




^ 5.4 X 10^ 



where the asymmetries a' for these modes are given in Table 2. 

Since each of these modes as well as several other possibilities have a 
different values of it is at least likely that a few instances of this kind of 
CP violation can be observed in the N'^^ ~ 10^ range. 

For comparison let us consider one of the CP eigenstate modes as in 
GLW, using the above numbers. In this case X = Xcp is a CP eigenstate. In 
particular, let us take the mode Xcp = Kstt^- Using the k = K* case we have 
as before a{K*) = 6.6 x 10"^ b{K*) = 10.0 x 10"^ In this case c(i^',7r°) = 
c{Ks'n-^) = 1.05 X 10"^. From this we get R = 0.24. Here 5 = so making 
a similar assumption as above concerning the phases, that sin ^ sin 7 ^1/2 
and cosCcos7 ^ we obtain a'{K*,KsTc^) ~ 0.12. Substituting this into 
equation (0) we obtain iV|r ^ 9.0 x 10''. 

In principle we can improve this by folding in information from all possible 
modes that are CP eigenstatesQ Some particular decay modes which are 
CP eigenstates are 



{3Ks, K^r], Ksp°, KsU, KsT]', Ksfo,Ks(f), 

Ksf2, vr+TT", K+K-} (53) 

Using the observed branching fraction in |Q these have a total branching frac- 
tion of about 5% (if the corresponding modes with are included, this total 
roughly doubles). Using the effective c obtained with the CP non-eigenstate 
method from all of these modes together, we get N^^{K*~ , Xcp) ~ 1.9 x 10^ 
which is roughly what we obtained in the single mode N^{K*~,K^p~) ^ 
2.5 X 10''. Thus, although the event rate for the various DCS modes is smaller 
than the CP eigenstate modes, the CP violating asymmetries are larger and 
iVg'^ is comparable. The latter case does have the following potential advan- 
tages 

^In particular we add the d(k,X) for CP even states to d{k,X) for CP odd states and 
vice versa 



20 



1. Each mode has a different strong phase difference and so it is perhaps 
more lilcely that one will have a total strong phase near ±7r/2 (i.e. 
maximal CP violation). 



2. Since a! tends to be larger, the bounds placed on 7 will tend to be more 
restrictive. 

In order to illustrate these two points, Fig. 2 shows the locus of allowed 
points on a plot of 7 versus h{K*') which would be obtained with the example 
above in the case of a K^'k" final state (sohd line) and a Ks'k^ final state 
(dashed line) assuming an ideal measurement. Clearly, in this instance the 
restriction on 6(A';) and 7 is much tighter for the K^t^~ final state than for 
the i^s7r° final state. 

5.3 Combining Information From Several Modes 

Let us now consider the degree of accuracy in the determination of 7 that 
may be obtained through the combination of various modes. To do this, we 
will do an illustrative calculation of the 90% confidence levels one can obtain 
on the 7 - h{k) plane with 7V|'' = 10^. 

To produce a scenario, we will consider some specific values of 7 and 
together with the branching ratios above and then calculate the values 
of d and d which would be relevant. If we now take these as being the 
experimental value we can then consider the likelihood as a function of 7 
and hiky 

The modes we will consider are (D^ K^ti~ , — > Kstt^, — > K~^p~ , 
K+a^, K^p^ and D° K*+n~). We will take Nb = 10^ and 

consider k = K* . It is important to note that in this example we consider 
statistical errors only. It is clear that in order to perform such a study, 
systematic errors must be under control. 

The exact experimental results will, of course depend on the strong phases 
and 7. For the purposes of illustration, let us assume that 7 = 60° and the 
strong phases: 



e(ir+7r) 



200 



30' 



10' 



1° 



1° 



e(ir*+7r-) 



20° 



40° 



50° 



(54) 
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Numerically, the values of d and d are given in Table 2 where the di, di 
entries are given in units of 10~^. 

In Fig 3a and 3b we show the likelihood contours as a function of 7 and 
b{K*). In Fig 3a we consider only the data from the i^+vr" and Kgn^ final 
states^. The solid curve shows the locus of solutions which explain the K~^7t~ 
data while the dashed curve shows the solutions which explain the Kgir'^ data. 
As can be seen there are four intersections which is the case in general when 
just two modes are considered. The contour regions correspond to 68% and 
90% confidence levels based on N^^ = 10^. In Fig. 3b we consider the 
situation where data from all six modes is used. In this figure, the solution 
for K~^7^~ is shown in solid, Ksn^ with short dashes, K~^p~with long dashes, 
K^Qi with the dash-dot line, Ksp'^ with the dash-dot-dot line and K*~^7t~ 
with the dash-dash-dot line. As can be seen from this figure, the correct 
value of 7 and b{k) is now identified with the only ambiguity being the four 
fold ambiguity in 7 between ±7 and vr ± 7. 

In Fig. 4 we have projected the likelihood from Fig. 3b onto the 7 axis 
where we have considered the case of 7 = 15°, 30°, 60° and 90° which are 
indicated by the curves peaked at those values of 7. In each of these cases, 
the 90% confidence interval is about ±10° about the solution. 

It should be realized that three body states K~^p~ , Kgp^ and K*^Tr~ can 
all lead to the common final state KsT^'^ir'. If one examines the distribution 
in phase space, then the vector resonances overlap to some extent and the 
channels will interfere with each other. In the following section, we will 
discuss how the additional information implicit in this situation can assist in 
extracting the value of 7. 

6 Using Three Body Decays 

Here we will consider the generalizations of the two approaches considered 
in section ^ to the case of a three body decay. First of all, we can consider 
the three body decay as consisting of a number of quasi two body channels 

•^Note that it has been suggested [p2[ that the combination of a DCS mode and a CP 
eigenstate mode (as in this example) offers some advantage because of the larger statistics 
that can be obtained by combining different eigenstate modes. This may be the case 
although how well a particular combination works also depends to a great degree on the 
strong phases. 
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which we can regard as distinct modes and find a solution for b{k) and 7. A 
second approach is to regard each point of the Dahtz plot as a distinct mode. 
We can then apply the inequalities eqns. ( PD| , |55D at each point. Since all of 
these inequalities must be true simultaneously, a very stringent bound can 
generally be placed on 7 and b{k). In fact we will argue that for at least some 
points this inequality is an equality so the limit given by such an argument 
should in fact give 7 and b{k). 

As an example we will consider in particular the case of D^, —>■ 
K^n^TT^. In this case the CBA decay D K^ti^'k^ has been experi- 
mentally studied by the E687 collaboration fl^. The data they obtain is fit 



to an amplitude to a general multi-channel 3-body decay form: 



M{W K+n-n^) = aoe'^" + ^ aiexp{i5i)B{a, b, c\r) (55) 

i 

where r is a label for the resonance and a,b and c are labels for the three 
final state particles which are permuted so that (a, b) forms the resonance 
that a given term represents. Thus, the function B is given by B{a, b,c\r) = 
BW{a, b\r)S{a, c) where if is the spin of r. 



BW{a,b\r) = -{sab - + iTr{Sab)mr) ^ 



1 if Jr = 

5(a,c) = { -2Pa-Pc if Jr = l (56) 

2(3(Pa-PcY -\Pa\''\Pcn if Jr=2 



where Sab = {Pa +Pfe)^ and Pa and Pc are the 3-momenta of a and c in the ab 
rest frame and 
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2 \ Jr+^ 



Tr{sab) = Trimi) | ^— — — — ^ — ^ \ (57) 

Here nir is the mass of the resonance and is the energy dependent width. 
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Thus, the general model used in |]T^ for a 3-body decay includes a number 
of resonance channels, each of which may have a different spin and phase 5i 
together with a non-resonant term (i.e. aoexp(i(5o)) 

Specifically, for —>■ K^tt^tt^ they use a non-resonant term together 
with three channels: 



M = aoexp{i6o) + a{K*^)exp{i6K'o)B{K+n-n°\K*°) 
+a{K*+)exp{i6K*+)B{K+7v\-\K*+) 
+a(p-) exp{iSp- )B{n-n^K+\p-) (58) 



In this analysis ||14|] the decay fractions^ for each channel and the relative 
phases were determined which we have summarized in Table 3. In our cal- 
culation we will use these decay fractions together with the global average 
for the total branching ratio, Br(^ K+tt^tt") = 13.9 ± 0.9% given in [§]. 
For the purposes of our illustrative calculations we will assume that these 
quantities take on their central values. 

The DCS decay — > K^7i~7t^ has not likewise been investigated exper- 
imentally, so for the purposes of making numerical estimates, we will use a 
model based on SU(3). First we write a decomposition of similar to that 

of eq. my. 



M{D'^ ^ ir+7r-7r°) = a'oexpiiS'o) + a'{K*'^) exp{iS'K,o)B{K+7T~7T'^\K*'^) 

+a'{K*+)exp{iS'K,+)B{K+7T%-\K*+) 
+a'{p-) exp{i6'p-)B{7i-7r^K+\p-) (59) 

Using the relations between the branching ratios to two body decays dis- 
cussed in section |5?1| , we can relate the magnitudes \a'{K*^)\, \a'{K*^)\ and 
\a'{p^)\ to \a{K*'^)\, |a(p~)| and |a(-ft'*+)| by regarding them as amplitudes 
for quasi- two body decays. In addition we set {a'^l = A^|ao|. Assuming exact 
SU(3) for the phases of the DCS channels, they will be: 

^ The decay fraction for a given channel X is the rate of D'~' K^tt^tt'^ if all 
channels aside from X are turned off divided by the total rate of — + K~^tt~tt^ 
with all channels present. Thus, for instance, the decay fraction through K*^Tr^ is 
Br{D K*^Ti^)Br{K*^ X+vr"). Due to interference effects, the decay fractions need 
not add up to 1. 
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^0 = '^0 = TT + 5k*o (5^.+ = 5p- 5'p- = 5k*+ (60) 

This model then predicts: 

Br{D^ K+n-TT^) = 7.8 x 10"^ (61) 

The only free parameters now are the strong phase difference ^ between 
k^D^ and B~ — >■ k~D , the overall branching ratios a{k) and b{k) 
and of course 7. If we use the estimates of a{K*) and b{K*) discussed in 
section ^, we can therefore obtain the Dalitz plots of the K~^7r~n^ final state 
which resulting from interference of D° and -D" channels. 

Fig. 5a shows such a plot with 7 = 90° and ^ = 70° which is intended 
to represent a sample of N'^ = 10^^. The upper plot represents B~ decays 
where the variables used are s = (ptt- + Pk+Y and t = (p^ro + Ppi-Y while 
the lower plot represents the CP conjugate. 

It is clear that in this case there is significant amount of CP violation. 
In Fig. 5b we have taken ^ = 160° and again it is clear that CP violation 
is present but here it is more subtle. There are about the same number of 
points in each of the plots but the distribution changes going from B~ to 
B+. 

In these plots, the three resonant bands are prominent, one could isolate 
the resonant contributions by fitting the plot to a model such as discussed 
above and thus obtain four modes (the three resonances plus the non-resonant 
contributions) to feed into the analysis of the previous section. Obviously 
there is more information implicit in these distributions. For instance, one 
would have additional constraints since one would know the relative strong 
phase difference between each of the D decay modes. 

In particular, for this model one has the magnitudes of four amplitudes 
each for B~ and B^ decay (giving 8 parameters) and three phase differences^ 
giving a total of 14 parameters. The three unknowns in this case are just 7, 
b{k) and ( so they are well over-determined. 

In order to get an idea of how many B mesons are required to extract 
useful information from this kind of data, let us consider what is required to 

^one cannot determine an overall phase so only the differences between the S's can be 
extracted from measurements 
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find a 3-a signal of CP violation for this system. In Fig. 6a we show a plot 
of the overall partial rate asymmetry for this mode for 7 = 90° as a function 
of the overall strong phase difference 77. In Fig. 6b we show with the solid 
line the number of (N^^) mesons required to give a statistical 3-a signal 
for the partial rate asymmetry. 

When the partial rate asymmetry is small, for instance in the case where 
1] = 160 in Fig. 5b, clearly iV^°" becomes large but of course, we are not using 



the information contained in the full distribution. Using the methods of JTS 
if we assume the CBA and DCS decays of the are understood we can 
construct an observable or system of weights of various regions of the Dalitz 
plot which is optimally sensitive to CP violation. Using this method, A'^^'^ is 
shown in Fig. 6b with the dashed curve; depending on 77, N'^^ varies between 
(0.35-2) X 10^ 

Another approach to obtaining 7 is to use the generalization of the bound 
given in eqn. (|30|). If one fits the Dalitz plots to continuous curves, then, 
regarding each point as a separate mode, we would have Qmin < Q as a 
function of the Dalitz plot variables. Since each lower bound must be valid, 
the maximum lower bound, Q = max{Qmin) will provide the most stringent 
lower bound on Q so therefore Q > Q. 

In fact it is not unreasonable to expect that Q = Q. To see this note that 
eqn. (p2D tells us that if Q = Qmin, then 



u = Zi + 2 cos^ 7 = M + 2-\/m COS 7 cos + 2 cos^ 7 (62) 
so therefore 



cos^ = — cos7/\/m (63) 

If eqn. ( |63|) is true anywhere on the Dalitz plot, then Q = Q and one would 
expect that this condition would apply on some curve on the dalitz plot. 
This is illustrated in Fig. 7a where the locus of points where Q = Qmin is 
shown for 7 = 60° with C = 0° (sohd line); C = 30° (dashed line); ( = 60° 
(dotted line) and ( = 90° (dot-dashed line). 

To extract Q using this method, it is useful to consider the function 

j^^^ _ fpalitz Q{Qmin{s, t) — q) ds dt ^^^^ 

Soaiitz ds dt 
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where the integral is over the allowed region of the Dalitz plot and 6{x) = 
if X < or 1 if X > 0. Thus f{q) is the fraction of the Dalitz plot such that 
Qmin > Q- For values of q Q from below, /(g) oc a/Q — q and so P{q) will 
linearly extrapolate its endpoint at q = Q. 

In Fig. 7b we show a plot of f{rQ) where r = 0.7 and 0.9 as a function 
of 7 for C = 0° (sohd); C = 30° (dashed); C = 60° (dotted) and C = 90° 
(dot-dashed). From this is apparent that for almost every combination of 
strong and weak phases roughly 20% of the dalitz plot has Qmin > 0.9Q. 

In Fig. 7c we show a graph of /^(g) for ( = 90° with Q = 0.25, Q = 0.5, 
0.75. As can be seen, in all cases a significant fraction of the Dalitz plot has 
Qmin close to Q and the curves extrapolate linearly to an endpoint at q = Q. 

Likewise we can extract b{k) by applying eqn. (^) to each point of the 
Dalitz plot. Here one has both an upper and a lower bound on b{k) so that 
if one considers the functions: 



g . (^h^ = I Dalitz ^ {bminjs, t) -b) ds dt 
S Dalitz dt 

= InaiuzOib-bmUs t)) ds dt 

^ ' iDaUtz dsdt ^ ^ 

the support of gmin must lie entirely below b{k) while the support of gmax 
must lie entirely above b{k). As with the bound in Q, the end point of the 
support in both cases is b{k). 

In Fig. 8 we show a graph of Qmin^b) and Qmax^b) as a function of b/b{k) 
in the case of —>■ K'^ir^n^. Here we take 7 = 60° and C = 90°. Again a 
significant portion of the Dalitz plot can be seen to give values of bmin and 
bmax close to the endpoint. 

If b{k) were determined in this way, the information could easily be com- 
bined with the branching ratio to a two body final state (e.g. K^tt^ or CP 
eigenstates) to obtain 7 as discussed in Case 1 of section |[ 

We can also apply the methods discussed in this section to the case where 
the decays to a 2-body final state (eg. D° — >■ K~^7r~ or Ksii^) but the 
parent decays to a 3-body final state, for instance D^Kstt^- In 

this case, of course the Dalitz plot variables are those of the parent B' decay 
but other than that, one could solve for 7 by fitting the the distribution to a 
series of amplitudes which could be treated as different modes or by finding 
Qmin as described above. 
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7 Summary and Conclusion 



In conclusion, we have considered a number of ways of observing the in- 
terference between ^ K~D^ and B~ K~D as a means of cleanly 
extracting 7. In general, this is achieved by observing final states common 
to and D decay. To fully determine 7 with such observations, several 
different decay modes must be observed or else the difficult to measure decay 
rate for K^D must be independently known. 

Thus if for one decay mode X the branching ratios for B~ K~X 
and B^ K^X are known, where the X results from the interference of 

and D channels and if one also knows B~ K~D then 7 can be 
determined up to an 8 fold ambiguity. The special case where X is a CP 
eigenstate was originally considered in 0. 

On the other hand if the branching ratio for B^ K^D^ is not known 
but there is a large degree of CP violation in the mode, then it may be 
possible to put a lower bound on sin^ 7. A particularly promising class of 
decays which could give large CP asymmetries are cases where D° — >^ X is a 
DCS decay. In such cases the enhancement of the decay X which is 

CBA is balanced by the enhanced production amplitude B^ —>■ K^D^ which 
is CLA. The two channels thus have roughly equal amplitudes therefore CP 
asymmetries can be large. 

If two or more modes Xi, . . . , X„ are measured, then we no longer need 
to know B~ — >■ K~D , it can be fit for along with 7. In the case of 2 
modes there is potentially a 16 fold ambiguity due to the need for solving 
of a quartic equation in sin^ 7. If 3 or more modes are considered, then 
there is only a 4-fold ambiguity since the system of equations for sin^ 7 is 
over determined. In our sample calculation it was found that for Nb = 10®, 
that the 90% bound on 7 found from combining several modes is roughly 15° 
which is typical although the actual precision depends on the values of the 
strong phase differences. 

This approach can be generalized to the case where the undergoes a 3- 
body (or indeed n-body) decay. If one considers each point in the dalitz plot 
to be a single mode, one can obtain a lower bound for sin^7 at each point. 
In general, one expects that the maximum lower bound is in fact equal to 
sin^ 7 so this method actually gives 7 up to a four fold ambiguity. The same 
method may also be applied to determine b{k) by obtaining upper and lower 
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bounds on b{k) at each point. Alternatively, fitting the Dalitz distributions 
to a resonant channel model provides enough information to obtain 7. 

Although throughout we have assumed that DD mixing is negligible, in 
the appendix we show how the effects of such mixing may be eliminated 
by using information about the time between the B~ and decays at the 
expense of increasing the statistical errors by about ^/2. 
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Appendix: The Implications of - Mixing 

In the discussion so far we have exphcitly assumed that D^D^ was negligible. 
In particular, since we often take advantage of interferences involving DCS 
decays which are 0(1%) of the interfering CBA decay, the total probability 
of mixing must be less than 0(1%) for the above formulation to remain valid. 
In this section we consider the generalization to the case where DD mixing 
may be present. 

We will argue that for final states which involve the DCS decay of the D^ 
(such as B~ K^[D^ K^ti^]) the effects of such mixing will be at most 
0(10%) on the rates (i.e. d{k,X)). 

In particular there are two possible ways to deal with mixing 

1. Using information on the the time between the B~ decay and the sub- 
sequent D^ decay, then the effects of possible mixing can be eliminated. 

2. If the parameters of DD mixing are known independently, then they 
can be taken into account in interpreting the time integrated data 

Indeed, if the mixing parameters and time dependent data is available, then 
one can in principle extract 7 from just one mode though most likely, the 
time dependence in the decay is too weak to make this a useful method. 
In order to examine the question, let us first define the standard parame- 



terization of this mixing as described in |T6[. We denote the mass eigenstates 
as 
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\Dl) = p\D'') + qW) \Dh) = - qm (66) 

where Dl and represent the hght and heavy eigenstates respectively. If 
we denote the ratio 



q/p = r,pe^"^ (67) 

where is a CP violating phase which could, in principle, be present in 
D° — D mixing and CP is also violated if Vgp ^ 1 though here we will 
suppose that would likely be the case for CP violation generated 

by physics beyond the SM. 

Let us denote by mn and rriL the mass of the Dh and Dl states respec- 
tively. Likewise, we denote by Th and the widths of these states. In 
terms of these quantities, we will use the definitions: 



m= , r = , // = m-zr/2. 

Am = niH - rriL, AF = F^f - F^, 

Am AF 
xd = -jr^ Vd = zd = xd + iyn = Ze'^, (68) 

and denote by \Dphys{t) > and \Dphys{t) > the time evolved state which is 
created at t = as D^f^y^ and D^^f^y^ respectively. These are thus related to 
the flavor eigenstates as: 



\Dphysit)> = g+(t)\D'' > +r,pe'''W >; 

\Dphys{t)> = 5+(*)l^>+r-V^^'^P°> (69) 

with g± are given by: 



g^{t) = e-^^*cosh(^zoT); g^{t) = e'"^' sinh{^ZDT). (70) 
with T = Vr)t. 
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A number of experiments have recently produced results which bound 



these mixing parameters. The E791 experiment at Fermilab obtains [[17 
Ar = 0.04 ± 0.14 ps'^ corresponding to y = 0.8 ± 2.9 ± 1.0% from the 
study of the life time ratio of —* K~tt~^ versus K^K~. CLEO has 

reportedm a 95% c.l. bound of \x'\ < 2.8% and -5.8% < y' < 1.0% where 



x' = xcos6k+tt- + y^^^^K+TT- and y' = —y sindK+w- + x cos6k+tt- obtained 
through the study of the time dependence of the decay D° — > K~^n~. Results 
from the FOCUS experiment ||19|| at Fermilab, again using the lifetime ratio, 



give y = 3.42 ± 1.39 ± 0.74% which interestingly is 2.5-cr from 0. 

In the Standard Model, -D mixing receives contributions from both 
short distance and long distance processes. The short distance mixing via 
the box diagram may be reliably estimated to be yn < = O(10~^). On 
the other hand, the long distance effects involve considerable uncertainty 
from the hadronic interactions involved. Calculations based on dispersion 
theory pO[ suggest that \zo\ < 10^^. On the other hand it has recently been 
suggested that life time differences driven by kaonic resonances in the 
vicinity of the mass could lead to y^ ~ 10~^. In the Standard Model, the 
CP violating phase will in all cases be negligible. 

It is also possible that — D mixing is driven by new physics beyond 
the standard model. Many examples of such mixing have been considered 
(see compilation in and in this case, additional CP violation through 
the phase angle is also possible. 

In the case we are interested in, B~ K~[D^ X], the initial decay 
is to a mixed state of — where both the and components have 
some amplitude to decay to X. Since X is small, it is justified to expand the 
mixing effects on the time dependent decay as follows: 



-^d(k,X) ^ (do(k,X) + di(k,X)T)e-^ 
dr 

4-d{k^X) ^ (do{k,X) + di{k,X)T)e-^ (71) 
dr 

Since do and do represent the decay rate at time interval, they will be 
identical with the expressions given in eq. (p2D. If we assume that there is 
no CP violation in the DD system, then di and di are given by: 
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^{di + di) = +2^ ac{X)Z ^^^Jac{X) sm(A + S) C0S7 + \/bc{X) sm(A - C) cos(7 + 0)| 
+2^bc{X)Z S^^Jbc(X) sin(A - 5) C0S7 + y'ac(X) sin(A + C) cos(7 + 0)| 

^(c?i-cli) = -2^Jac(X)Z !^^Jac{X) cos(A + 5) sin7 + \/bc(X) cos(A - C) sin(7 + 0)| 
+2^6c(X)Z |^6c(X) cos(A - 5) sin7 + ^ac(X) cos(A + C) sin(7 + (/))| 



In order to estimate the relative effect of the mixing, let us assume that 
Z ^ 0.01 and recall from our rough estimates that 

^«^-0.01 (73) 

a c(a) 



and thus from eq. (^) we can estimate: 



^ ^ i ^ 0.1 (74) 

do do 

Note that in Eqn. ([72|) the time dependent term depends on both ( and 
6 rather than just the sum ^ and, in addition, it depends on the mixing 
parameters. The simplest way to extract 7 given time dependent data is 
therefore to try and obtain the values of do and do which do not have these 
extra dependences reducing the problem to the same situation as if mixing 
were not present. 

This can be accomplished through weighting the data with 



(72) 



wo{t) = 2- t 



(75) 



so that 



do = / [d{T)wo{T)]dT] do = / [d{r)wo{T)]dT 
Jo Jo 



(76) 
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Using this method more data would be required to obtain the same sta- 
tistical results as in the unmixed case. In the unmixed case where di = one 
could obtain do more effectively by taking the time integrated rate. Thus in 
the unmixed case if a measurement of dQ is based on n events, the uncertainty 
in do is given by: 

^ = — {no mixing) (77) 

uq n 

In the mixed case, using eqn. (|T6D the uncertainty is 



{Mo)' 
dl 




{with mixing) (78) 



From Eqn. (0) this means that roughly twice the data is needed to have 
the same statistical power as in the unmixed case. In order to gauge the 
precision of time measurement required, we can smear out the distribution 
in eq. (|76D with a Gaussian resolution function of the form 



r(r,r') oce 5^ (79) 

where r is the actual time of the decay, t' is the measured time of the decay 
and a is the resolution (all in units of 1/To)- Since r is symmetric under 
r ^ r', the fact that w is linear in r implies eq. ([76| ) will still be true for r' 
but now the error is: 



(Ado) 



d\ 



2 




(80) 



As can be seen, the number of events required is not adversely effected if 
c" < l/^D but will be significantly degraded otherwise. 

In the above strategy, it was assumed that the mixing was relatively 
unconstrained except for the assumption that di is smaller than do which 
seems justified. If the rate of DD mixing is near the current experimental 
bounds, however, then there is another possible way to obtain 7. If one 
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measures {do, di, do, di} and knows all of the mixing parameters, then there 
are just four unknown parameters that need to be fitted for: {7, 6, C, and 
b}. We therefore have four equations in four unknowns and can solve for 7. 
To obtain di we can use the weight function: 



1 



(81) 



It is also possible to learn about 7 using time integrated data. This has 
been previously considered in 



23| . The corrections due to mixing in the 



approximation we are using are: 



d = do + di 

d = do + di (82) 

In the case that the DD mixing parameters are already determined, if n 
modes are measured, there are 3 + n unknowns: (, 7, b and for each mode 
Xi, 6i. Each mode provides two pieces of information, d and d, leading to a 
total of 2n measurements. To solve for the unknowns at least 3 modes are 
thus required (six equations in six unknowns). If the mixing parameters are 
unconstrained, this adds four more unknowns: {x, y, Vgp, (f)} and in principle 
at least seven modes would be required. 

On the other hand, using such time dependent information is probably not 
the best way to measure DD oscillation parameters. Perhaps the best way 
to proceed is to consider di and di to be corrections to the time integrated 
data. From eqn. (|72D one can take the bound on x and y together with an 
estimate of b and bound these terms. As indicated, this gives about a 10% 
correction which, as discussed in leads to about a 15° error in 7. This 
is similar to the errors typically obtained in the unmixed case for Nb = 10^ 
and so to gain improvements in the determination 7 with larger numbers 
of Nb in this scenario one would both have to improve the bounds on DD 
mixing (or indeed discover it). Of course since measurements of DD mixing 
can often be made at a 5 facility, it would be natural for such improvements 
to occur at the same time. 



References 



34 



[1] L. L. Chau and W.-Y. Keung, Phys. Rev. Lett. 53, 1802 (1984). 

[2] M. Gronau and D. Wyler, Phys. Lett. B265 (1991); M. Gronau and 
D. London., Phys. Lett. B253, 483 (1991). 

[3] L Dunietz, Phys. Lett. B270, 75 (1991). 

[4] L Dunietz, Z. Phys. C56, 129 (1992). 

[5] L Dunietz, "CP Violation with Additional B Decays", published in B 
Decays, S. Stone ed. (World Scientific, Singapore, 1992). 

[6] D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 78, 3257 (1997). 

[7] T. E. Browder, K. Honscheid and D. Pedrini, Ann. Rev. Nucl. Part. Sci. 
46, 395 (1996); see also M. Neubert, V. Riekert, Q. P. Xu and B. Stech in 
Heavy Flavors, A. J. Buras and H. hinder (World Scientific, Singapore, 
1992). 

Particle Data Group, Z. Phys. C3, 1 (1998). 

[9] A. Ah and D. London, Eur. Phys. J. C9, 687 (1999); 

[10] N. Sinha and R. Sinha, Phys. Rev. Lett. 80, 3706 (1998). 

[11] H. Yamamoto, Proceedings, Fourth International Workshop on Particle 
Physics Phenomenology, Kiohsiung, Taiwan (1998) [hep-ph / 9812279 . 



[12] A. Soffer, Phys. Rev. D60, 054032 (1999). 
[13] A. Soffer, |hep-ex/98010T8 . 



[14] P. L. Frabetti et al, Phys. Lett. B331, 217 (1994). 

[15] D. Atwood and A. Soni, Phys. Rev. D45, 2405 (1992); J. F. Gunion, 
B. Grzadkowski and X. He, Phys. Rev. Lett. 77, 5172 (1996). 

[16] Z. Xing, Phys. Rev. D55, 196 (1997); T. Liu, Tau Charm Factory Work- 
shop, Argonne IL (1995). 

[17] E. M. Aitala et al [E791 Collaboration], Phys. Rev. Lett. 83, 32 (1999). 



35 



[18] R. Godang et al. [CLEO Collaboration], |hep-ex/UUUlU60 



[19] J. M. Link et al. [FOCUS Collaboration], |hep-ex/0004034 



J. F. Donoghue, E. Golowich, B. R. Holstein and J. Trampetic, Phys. 
Rev. D33, 179 (1986) 

[21] M. Gronau, Phys. Rev. Lett. 83, 4005 (1999). 

[22] H. N. Nelson, |hep-ex/9908021| . 

J. P. Silva and A. Soffer, Phys. Rev. D61, 112001 (2000). 



36 



Table 1: Cabibbo allowed (CBA) and doubly Cabibbo suppressed (DCS) 
modes of D^{D^). BR's for CBA modes are taken from PDB |^ while 
those for DCS are predicted, using the model given in the text, except for 
D K^ir^ wherein the measured value is shown. 



Mode 


Br{D^ final state) 


Br{D'^ final state) 


K+TT- 


(2.9 ± 1.4) X 10^^ 


3.83 X 10-2 


K+p- 


3.8 X 10-^ 


10.8 X 10-2 




7.0 X 10-5 


7.3 X 10-2 




8.3 X 10-^ 


5.0 X 10-2 
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Table 2: The branching ratios for the combined decay —>■ K*^D^ fol- 
lowed by the decay of to the modes given below using the parameters 
considered in Section 6 with 7 = 60° and the given strong phases ^j. di and 
d are given in units of 10"*^ and a' is the partial rate asymmetry. 



Mode 


di 


di 


\{di + di) 


a' 






91 


75 


83 


0.096 


10 




842 


740 


791 


0.064 


20 


K+p- 


289 


159 


224 


0.288 


30 


K+a^ 


203 


90 


146 


0.383 


40 


Ksp' 


333 


391 


362 


0.081 


200 




97 


34 


65 


0.477 


50 



Table 3: 
in 111. 



The parameters of the model for —>■ K tt+tt*^ decay obtained 



channel 


decay fraction (%) 


6r - 5p (deg.) 


non-resonant 
K*^ 
K*+ 


10.1 ±3.3 ±3.0 ±2.7 
16.5 ±3.1 ±1.1 ±1.1 
14.8 ±2.8 ±4.9 ±0.3 
76.5 ±4.1 ±2.2 ±4.9 


-122 ± 10 ±21 ±2 
-2 ±12 ±23 ±2 
162±10±7±4 
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Figure Captions 



Figure 1: Each of the sohd hnes shows the locus of points in 7 versus Ui of 
allowed solutions given Zi = 1.5 for i/i = (outer curve), 1 (intermediate 
curve) and 2 (inner curve). The boxes indicate the inequalities 
Eqns. (|30|J3|). 

Figure 2: The solid line shows the allowed points on the 7 — h{K*) plot 
given a measurement of B~ K*D^ followed by K^tt^ (and charge 

conjugate) using the estimated values of Section 5 eqns. ( ^ , ^ ) assuming 
7 = 90°. The dashed line shows the results which would follow in the same 
case given a measurement of —>■ KgH^ in the final state. 

Figure 3: (a) The likelihood distribution is shown as a function of 7 and 
b{K*) assuming that N^'^ = 10^ with the branching ratios considered in 
Table 2 and assuming only the K^tt" and KsTt^ modes are measured. The 
outer edge of the the shaded regions correspond to 90% confidence while 
the inner edge corresponds to 68% confidence. The solid lines show the 
locus of points which give exactly the K~^7r~ results while the short dashed 
curve shows the points which give the KsTT^ results, 
(b) The likelihood distribution as in Fig. 3a is shown assuming all of the 
modes in Table 2 are used. The solution for the K^tt~ data is shown with 
the solid curve; that for the Kst^^ data is shown with the short dashed 
curve; the one for the K^p~ data is shown with the long dashed curve; the 
one for the K^ai data is shown with the dash-dot curve; the one for the 
Ksp^ data is shown with the dash-dot-dot curve and the solution for the 
K*~^'n'~ data is shown with the dash-dash-dot curve. 

Figure 4: The ratio between the the likelihood distribution and the 
maximum likelihood is shown as a function of 7 with the parameters as in 
Fig. 3b except 7 is taken to be 15° (dashed curve); 30° (solid curve); 60° 
(dotted curve); 90° (dash-dot curve). 

Figure 5: (a) Dalitz plots for the decay products of the in the decay 
chain K*^D^ followed by K'^tt~tt^ (upper plot) and the 

charge conjugate process (lower plot). We use the model described in 
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section |^ where we take 7 = 90° and ^ = 90°. The plots represent the 
results given Nb = 10^°. 

(b) Dalitz plot for the same system as in Fig. 5(a) with ^ = 0°. 

Figure 6: (a) A plot of the partial rate asymmetry for the system in Fig. 5 
as a function of 7. 

(b) The Nb required for a 3 — a signal as a function of ^ using only the 
PRA (solid line) and using the optimal observable (dashed line). 

Figure 7: (a) The locus of points on a dalitz plot where Qmin = Q for 
7 = 60° and C = 0° (solid line); ( = 30° (dashed line); C = 60° (dotted line) 
and C = 90° (dot-dashed line). 

(b) A plot of f{rQ) as a function of 7 for ^ = 0° (sohd curves); for ^ = 30° 
(dashed curves); for C, = 60° (dotted curves) and for ( = 90° (dash-dot 
curves). In each case the lower curve corresponds to r = 0.9 and the upper 
curve to r = 0.7. 

(c) A plot of P{q) is shown as a function of g for Q = 3/4 (solid line); 
Q = 1/2 (dashed line) and (5 = 1/4 (dotted line) where C = 90°. 

Figure 8: A plot of g^irX^) (solid line) and g^axi^) (dashed line) line is 
shown as a function of b/b{k). Here we have taken Q = 3/4 and C = 7r/2. 
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Figure 1 
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Figure 3a 
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Figure 3b 




15 30 45 60 75 90 
gamma deg 



44 



Figure 4 
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Figure 5b 



B- gamma= 90.0 zeta= 160.0 
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Figure 6a 
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Figure 7a 
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